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The information in quantum computers is often stored in identical two-level systems (spins or 
pseudo-spins) that are separated by a distance shorter than the characteristic wavelength of a 
reservoir which is responsible for decoherence. In such a case, the collective spin-reservoir interaction, 
rather than an individual spin-reservoir interaction, may determine the decoherence characteristics. 
We use computational basis states, symmetrized spin states and spin coherent states to study 
collective decoherence in the implementation of various quantum algorithms. A simple method 
of implementing quantum algorithms using stable subradiant states and avoiding unstable Dicke's 
superradiant states and Schrodinger's cat states is proposed. 

PACS numbers: 03.67Lx, 03.65Yz, 03.67.-a 



I. INTRODUCTION 



In future quantum information processing systems, 
quantum bits (qubits) of information are likely to be 
stored in identical two-level systems (spins or equivalent 
spins) that are spatially separated to simultaneously al- 
low independent control of single qubits and mutual in- 
teraction between qubits. If a system of n spins interacts 
with a reservoir whose characteristic wavelength is longer 
than the inter-spin distance, the system interacts collec- 
tively with the bath. For instance, electron-spin-based 
quantum computers 1] have an inter-spin distance on the 
order of 10 to 100 nm, which can be far shorter than the 
wavelength of a coupled low temperature phonon reser- 
voir. Another example is a nuclear-spin-based quantum 
computer using doped impurities 2] or a crystal lattice 
[3, y in which a remote paramagnetic impurity interacts 
with a nuclear spin system by a long-range dipolar cou- 
pling . If the distance between qubits is much smaller 
than their separation from the impurity, the interaction 
can also be treated as collective. 

In theory, if a system-bath interaction exhibits a high 
degree of symmetry, computation can be performed in a 
decoherence-free subspace (DFS) that leaves the logical 
qubits immune to the influence of the bath 0, Q, H, 0, 
liol ITU IT^ ITsf . However, if one restricts logical states 
to this subspace, it is difficult to design simple gates to 
implement quantum algorithms. Consider, by compari- 
son, a standard model of quantum computation, which 
labels logical states (computational basis states) by the 
tensor-products of the two-level qubit states, and imple- 
ments the unitary evolution for a given algorithm by a 
sequence of controlled-NOT (CNOT) gates and single- 
qubit rotations. These fundamental gates are physically 
convenient, as they can be created simply through pair- 
wise interactions between qubits and local control pulses 
(i.e., by exploiting the innate tensor product structure 



of the Hilbert space) ^. By contrast, if logical states are 
defined within a DFS, the control sequences are compli- 
cated by the requirement that the state does not wander 
out of the DFS during evolution. Thus, in the interest of 
simple gate construction, we examine decoherence in the 
context of the standard computational basis and the full 
Hilbert space. 

Consider the unitary evolution of a quantum computer 
during the course of a given algorithm. Ideally, the state 
is described by a trajectory in the 2"-dimensional com- 
plex Hilbert space. If a collective system-environment in- 
teraction is present, the instantaneous decoherence rate 
can be parameterized by the "proximity" of the ideal tra- 
jectory to the decoherence-free subspace, allowing us to 
infer how the decoherence rate evolves over the course 
of the algorithm. Generally, the end result of a quan- 
tum algorithm is encoded in the measurement of the final 
state; there is thus some freedom to choose an initial state 
and/or trajectory that minimizes the average proximity 
to the DFS. In this fashion, one can attempt to minimize 
the impact of decoherence due to the collective system- 
reservoir interaction while maintaining the simplicity of 
the tensor-product structure of the full Hilbert space. As 
a compromise between a naive application of the stan- 
dard model and logical encoding in a DFS, this approach 
might be useful to implement algorithms on a prototype 
quantum computer, where the number of qubits is too 
modest to implement quantum error correction. 

In this paper, we introduce physical models for collec- 
tive spin-reservoir interactions, and examine relevant pa- 
rameterizations of the decoherence rate during the time 
evolution of specific useful quantum algorithms. It is 
shown that judicious choice of the initial state or gate se- 
quence can often reduce the decoherence rate parameter 



^ Note that we simplify matters by neglecting fault-tolerant imple- 
mentation and system-environment interactions that break per- 
mutation symmetry. 
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significantly. 

A collective spin- reservoir interaction that is responsi- 
ble for a longitudinal relaxation process (Ti process) can 
be modeled by 

Hi = Y,fih[alj- +aj+y (1) 

i 

Here, cii (a|) is an annihilation (creation) operator of 
the bosonic reservoir and J+ (J^) is a collective raising 
(lowering) operator of the spin system: 

1=1 1=1 

where a'^ = a1±. icr^ are Pauli operators for spin £. 

In this model, the longitudinal relaxation time T\ can 
be much shorter than that for a single isolated spin; this 
phenomenon is known as superradiant decay in atomic 
physics fT3|. In contrast, the so-called subradiant states 
have negligible 1/Ti. To characterize this longitudinal 
relaxation rate, it will be useful to define simultaneous 
eigenstates of the squared total angular momentum 
and the z-component of the total angular momen- 
tum J, notated symmetrized spin states, and discussed 
further in Sec. Ill Bl By projecting the ideal trajectory 
onto the symmetrized states, we can calculate a time- 
dependent parameter to estimate the T\ relaxation rate. 

The collective spin-reservoir interaction that we adopt 
to model transverse relaxation i^T-i) processes is 

ni^^hg,alaj\ (3) 

i 

Here is the z-component of the total angular momen- 
tum J, 

N n 

^^=E^-2^-| ■ (4) 

In order to gain a physical picture for the collective 
Ti and T2 decoherence behaviors in an N spin system, 
it is convenient to use a spin coherent state which is a 
linear superposition of symmetrized spin states with the 
same total ang ular momentum j and the third quantum 
number a |l5i The Q{a) representation of the spin 
state in terms of spin coherent states can easily identify 
a macroscopically separated linear superposition state 
which is subject to strong collective decoherence. 

We will examine for modest n the collective deco- 
herence properties of two representative quantum algo- 
rithms: the Deutsch-Josza algorithm [l?.] and quantum 
search . We find a few cases in which the system state 
must pass through vulnerable states, such as Dicke's su- 
perradiant states and macroscopically separated linear 
superposition states, and so is vulnerable to rapid longi- 
tudinal relaxation and transverse relaxation due to large 
instantaneous coupling to the external reservoir. We 
propose that simple modifications to the algorithms can 
avoid this problem. 



II. COMPUTATIONAL BASIS STATES, 
SYMMETRIZED SPIN STATES AND SPIN 
COHERENT STATES 

In this section, we introduce three sets of basis states. 
The first set is a standard computational basis formed 
by the tensor product of n two-level systems, and is 
most natural in discussing gate-by-gate implementations 
of various algorithms. The second basis set consists of 
symmetrized states that are convenient to calculate col- 
lective longitudinal and transverse relaxation rates. The 
third is a set of spin coherent states that is suitable to 
discuss Q{0, ip) representation. 

A. Computational basis states 

In accordance with the standard model describing 
the state of n spin-1/2 qubits, computational ba- 
sis states are defined as simultaneous eigenstates of 

f t2 t2 t2. jz jz jz\. 

W 1 J ^^2 I ■ ■ ■ ' '^n I 1 I "^2 J ■ ■ ■ ; '^n I ■ 

•^|ji,i2,-- - ,j„;mi,m2, • • • , m„) 

= + l)lii, j2, ■ • ■ ,jn;mi,m2, ■ ■ ■ ,mn) (5a) 

J^\jlJ2,■■■ ,j„;mi,TO2, • • • ,TO„) 

= j2, ■ ■ ■ , j„;mi,TO2, • • • ,m„) , (5b) 

where ji — 1/2 and to^ = ±1/2. As all the ji values are 
fixed for a collection of two-level spins, we omit them in 
denoting a computational basis state. Thus, the Hilbert 
space is spanned by the 2" basis states |mi, TO2, . . . , mn). 

B. Symmetrized spin states 

1. Defimtion 

Symmetrized spin states are defined as simultaneous 
eigenstates of ^J^ , J2 , ■ ■ ■ , J^', , J^^- 

Ji Iii,j2, • • . ,j„;j,m) = + 1) j2, ■ • ■,jn;j,m) 

(6a) 

J^\jl,j2,---, jn ; j, m) = j (j + 1) I Jl ) j2 , • ■ ■ , jn \j,m) 

(6b) 

J" lii, J2, . . . , j„; j, m) = m | ji, j2, • • ■ , j™; j, m) (6c) 

We again omit the ji labels from the symmetrized 
states. Note that j € {0, 1, 2, . . . , n/2} for n even 
and j E {1/2, 3/2, . . . , n/2} for n odd, and m E 
{—j, — j + 1, . . . , j}. Generally, each symmetrized state 
|j, TO, a) is d-fold degenerate; a E {1, 2, . . . , d} encapsu- 
lates all other quantum numbers required to distinguish 
each basis state. Symmetrized spin states can be lin- 
early expanded in terms of computational basis states via 
Clebsch-Gordan coefficients, whose properties are sum- 
marized in Appendix El 



3 



2. Collective longitudinal relaxation rate 

To gauge the impact of the longitudinal relaxation pro- 
cess described by Eq. one could solve the Heisen- 
berg equation of motion for the joint density matrix of 
the system and reservoir, under the influence of the sys- 
tem and bath internal Hamiltonians, any externally ap- 
plied control fields, and the system-reservoir coupling. 
However, such an analysis is prohibitive due to the large 
number of degrees of freedom. To simplify, we view the 
system-reservoir coupling as a perturbation that causes 
transitions between system spin states. In the absence 
of any relaxation process (i.e., to zeroth order), the sys- 
tem evolves in a unitary fashion according to the algo- 
rithm of interest. The transition rate 7 due to system- 
environment coupling is calculated in Appendix^ this 
rate varies throughout the time evolution, as it is depen- 
dent on the system state. The transition rate serves as a 
parameter to characterize the instantaneous longitudinal 
relaxation rate. 

For the system-bath interaction described by Eq. Q, 
it is convenient to expand the system state in the sym- 
metrized basis, as the environment acts on the spins 
through the collective raising and lowering operators J* . 
The transition matrix elements for these operators are 



(j,m± l,a| |j, m, a) = {j T m) {j ± m + 1) . (7) 



It is assumed that the bath is cold and resides in the 
ground state; the downward transition is thus relevant. 
The factor in Eq. ((TJ represents the amplification in the 
collective relaxation rate relative to the transition rate 
7o for a single spin. Note that for n even, the states 
\j = n/2,m = 0) and \j = n/2,TO = 1) have the highest 
value of 7: 



7 = + m){j - TO + 1)70 



n f n + 1 



70 



(8) 



The longitudinal relaxation rate per spin (that is, 7/f^) 
is thus enhanced by a factor of approximately n/4 com- 
pared to that of an isolated single spin. For a quan- 
tum register of 4000 spins, the longitudinal relaxation 
rate would be increased by three orders of magnitude. 
This enhanced relaxation rate is known as superradi- 
ant decay 0|. On the other hand, the states \j,—j,a) 
(j — n/2—1, . . .) have zero longitudinal downward relax- 
ation rate. If a reservoir cannot provide enough energy 
to the spin system to cause an upward transition, those 
states are completely stable against a process and 
are thus referred to as "subradiant states" . Eq. iQ also 
shows that the collective longitudinal relaxation rate de- 
creases with the total angular momentum j. Since the 
downward transition matrix element for to = scales 
j{j + 1), the state |j = 1,to = 0) has a relaxation 
rate more than six orders of magnitude smaller than the 
superradiant state \j = n/2 = 2000, to = 0). 



An arbitrary spin state can be expanded as a lin- 
ear combination of an orthonormal complete set of sym- 
metrized states, 



^ ^ Cj.m,a\ii 



TO, a) 



(9) 



The probability density of eigenvalues j and to for a given 
state are 

a 

= ^\{j,m,a\ip)f 

a 

= {fp\^\j,m,a){j,m,a\ilj) . (10) 



Since the longitudinal relaxation rate is dependent on 
only j and to but is independent of a, the probability 
of finding a given state \ip) in the subspace designated 
by (j, TTi) is sufficient for our purpose of studying the 
longitudinal relaxation. 

The overall normalized longitudinal relaxation rate for 
an arbitrary state \^) is thus given by 



J_ 

70 



= ^^P(j, TO)b-(j + l)-TO(TO-l)] . (11) 



3. Collective transverse relaxation rate 

The collective transverse relaxation rate is calculated 
by a model for a quantum nondemolition (QND) mea- 
surement of , which is summarized in Appendix [CI 
The off-diagonal element \j,m,a){j' ,m' , a'\ in the den- 
sity matrix decays with a rate proportional to |to — to'| 
by such a QND measurement process. 

Since the transverse relaxation rate is only dependent 
on the quantum number to, we can expand an arbitrary 
71-spin state \tp) in terms of the computational basis states 

IV') = X! X!' ' ' X! (^mim2-m„ \mim2- ' ' TO„) . (12) 

After the T2 relaxation process, the (mixed) state density 
operator is 



mi,--- ."III m', .m' 



p = y . E 
i,--- 

xe-^°l™-"'lVi,-- - ,to„)(to'i,... ,to:,| (13) 

where to = ^Yi!i=\ '^ii = Sr=i ^^"^ -^0 is a Ti 
relaxation rate for a single spin. 

The fidelity of the mixed state p compared to the initial 
pure state j-)/') is adopted as a measure for the degradation 
of the state : 



(14) 
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If we substitute Eqs. (221 and into Eq. we 
obtain 




(15) 



where Fij = /o|r7ii — mj\. We approximate e ~ 1 — 
Fijt for a short evolution time {t <C If we use the 

relation X]i=o ^ S?=o^ I'^iPl'^j P — which is obtained 



from Tr(p) = 1, we have 

2"-l 2"-l 

2 



2"-l 2"-l 

^ = 1- 9 E E |c,nc,f/^o|m.-m,|t 



i=0 j=0 



1 - Irt. 

2 



(16) 



The normalized T2 relaxation rate for an n-spin system 
is then given by 



^ 2"-12"-l 



(17) 



Note that the states \j — 0,m = 0,0;), with a d — 
"-'/[(§ + !)■§'] fold degeneracy, form a decoherence free 
subspace (DFS). 



C. Spin-coherent states 

1. Definition 

Spin-coherent states \9, ip) are defined as eigenstates of 
the following non-Hermitian operator |19| |: 



J^sin6' + e"^cos2 [ -] J- 



e'^sin(^-j J+J =0. 

(18) 

Eq. H18() resembles the definition of a coherent state of 
harmonic oscillator, {a — a)\a) = 0, where a is a complex 
eigenvalue. Spin-coherent states in the j — n/2 subspace 
can be mathematically constructed by the rotation opera- 
tor Rg ^ acting on the ground state \j — n/2, m = —n/2) 
El: 



\e,v) = Re^^ -,--) (19) 
Rf,^^ = e»e(J^ Bin cosy) (20) 

As Re^ip commutes with J^, \9,Lp) is an eigenstate of 
with the same eigenvalue j(j -I- 1) as |n/2, —n/2). There- 
fore, |0, (/?) can be expanded as a linear combination of 
the completely symmetric states |j = 71/ 2, m). 



1/2 



rn— — nl2 



(l-H |t|2)"/2 \m + 



2'"^)' (21) 



where r 



^v'tan(l) . 



Spin-coherent states are not mutually orthogonal. The 
inner product of two spin-coherent states is given by 



|(0',^'|0,^)|2=cos2"($) 



(22) 



where $ is the angle between \9,if) and and is 

given by 

$ = cos~^ [cos 6 cos 9' + sin 9 sin 9' cos{ip - (p')] . (23) 
Spin coherent states satisfy a completeness relation: 

/. (24) 



(n+1) / ^\9,^){9,^\^Y.\^,m){^,m 



As the coherent states span the j = n/2 subspace and 
are linearly dependent, they form an over-complete set. 

According to Eq. ((T^ . spin coherent states are pro- 
duced by rotation of the ground state |n/2, —n/2), which 
is a product state without any entanglement. The oper- 
ator in Eq. (|20|l can be broken into a product of rotation 
operators for each individual spin: 



^iO^Ji sin ^ — J^ cos I/?) ^ fJ-^i'J^ ^~'^2 ^) 



pid{J^ sin :p~J^ cos ip) 



(25) 



The spin coherent states generated by this rotation op- 
eration are also unentangled product states. 



2. Q representation 

We introduce the Q-representation, defined as the 
overlap between a given state \tjj) and a spin-coherent 
state 16*, 95): 



Q{9,^) = \{9,^m\ 



(26) 



4, m = 


4)- 


Hj 


= 4, m = 


-4)), 


A,m — 


3)- 


Hj 


= 4, m = 


-3)), 


4, m = 


2)- 


Hj 


= 4, m = 


-2)), 


4, m = 


1>- 


Hj 


= 4, m = 


— 1)), and 



FigureOshows the Q{9, ip) representation of the following 
states: 

|j = 4,m = 0). 

We see a non-classical oscillation in the x-y plane due 
to interference between the superposed symmetric states. 
As the difference between m values of the components in- 
creases, the number of oscillations increases as well. As 
discussed in Sec. Ill B 31 the states with rapid oscillations, 
and thus widely separated m values, are particularly sus- 
pectable to T2 relaxation. In this fashion, the structure 
of the Q-function provides a measure of the collective 
transverse relaxation rate. 
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3. Spin coherent states for all j manifolds 

So far, we have discussed spin coherent states for the 
completely symmetric subspace {j = n/2) through Eq. 
(|19|l . We can define spin coherent states for arbitrary 
values of (0 < j < n/2 — 1) by applying the rotation op- 
erator Rs^ip to the subradiant states \ j < n/2, m = — j, a) 
[2(tI |. In this way, the Hilbert space can be partitioned 
into separate spheres with different radii j, as shown in 
Fig. El 

III. APPLICATION TO SELECTED QUANTUM 
ALGORITHMS 

In the previous section, we introduced three basis sets 
to aid in characterizing the decoherence during the execu- 
tion of representative quantum algorithms. Instances of 
two quantum al gori thms will be analyzed: the Deutsch- 
Jozsa algorithm [13 and Grover's data search algorithm 

M 

A. Deutsch-Josza algorithm 

1. Standard implementation 

The Deutsch-Jozsa algorithm consists of the following 
four steps as shown in Fig. 0] 

Stepl : Preparation of an n-qubit control register and 
a one-qubit work register. 

1^1) = |O)c0|l)^ 
Step2 : Hadamard transformation. 

1 ^""^ 1 

^ x—0 ^ 

Step3 : /-controUcd-NOT gate operation. 

1 ^""^ 1 

^ x—Q ^ 

Step4 : Hadamard transformation. 

y=0 X=0 



Xn-1, x„-2, ■ ■ ■ ,Xo be the binary digits of x. Let Xi rep- 
resent the logical complement of xf, the symbols V and 
A are logical OR and AND, respectively. The function 
f{x) can be expressed as 

f{x) — {xi A X2 A 2:3 A X4 A X5 A A xy A Xs) 
V {xi A X2 A X3 A X4 A X5 A xe A XT A xs) V . . . (27) 

Clearly, there will be 128 bracketed expressions for any 
balanced function with n = 8. Each expression can be 
computed in sequence in a quantum circuit by flipping 
the complemented bits^, performing a CNOT gate with 
all n qubits as control bits and the work qubit as the 
target, and again flipping the complemented bits. Here, 
we consider each of these 128 blocks as a single compu- 
tational step; along with the Hadamard transformations, 
this leads to a total of 131 steps in the algorithm, includ- 
ing the initial state^. 

For the f{x) described above, we have calculated the 
projection of the state after each of these 131 steps onto 
the symmetric state basis; here, we only show in Fig. |31 
the projections at the selected steps indicated in Fig. 21 
The normalized relaxation rates can then be calculated 
at these instants from the projections. In Fig. O the 
projections onto the spin-coherent states is shown. The 
initial state j'i/'i) has P{j, m) — dj_i6m,~4: , which is trans- 
formed to a spin coherent state |-02) — \0 — ■§ , = O) hy 
Hadamard transformation. This spin coherent state has 
nonzero P{j,m) values only in the j — n/2 subspace 
and is peaked at m = 0, so that the state is vulnerable 
to superradiant decay. After implementation of half of 
the /-controUed-NOT gates (|^3_2)), the Q-function ex- 
hibits a macroscopically-separated linear superposition, 
as shown in Fig. The state \tp3-3) immediately be- 
fore the second Hadamard transformation is again a spin 
coherent state with an opposite phase \9 = 7r/2, ip = t:) 
compared to the state |V'2)- 

Figure [T^V shows the evolution of the normalized Ti 
relaxation rate at each of the 131 computational steps in 
the implementation of this algorithm. The Ti relaxation 
rate is highest immediately after the first Hadamard 
gate and immediately before the second Hadamard gate, 
where the corresponding states |-!/'2) and |'03-3) are 
spin-coherent states with a superradiant relaxation rate 
7/70 = 18. This relaxation rate should be compared 
to the overall Ti relaxation rate 7/70 = 4 for the case 
that n = 8 spins are in the same spin-coherent state but 
interact individually with an external reservoir. 

The normalized T2 relaxation rate in the implemen- 
tation of this algorithm jumps from zero for the ground 



The /-controUed-NOT gate operation could be per- 
formed by expressing /(x) in disjunctive-normal form 
and adding CNOT gates with multiple control bits for 
each expression. For example, consider a balanced 
function /(x) that acts on n = 8 qubits, yielding 
one if the parity of x is odd and zero otherwise; let 



^ The reader may note that a much simpler, computationally ef- 
ficient circuit exists for the example provided. However, the de- 
composition described above is perhaps more general, and it is 
applicable to any balanced function. 

^ The consequences of this arbitrary partitioning of the algorithm 
evolution are discussed in Sec. IIVI 
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state to 7/70 — 1-51 for the spin-coherent state, and stays 
at the same rate for the remaining steps; this is due to 
the fact that the 128 n-control-bit CNOT gates used to 
implement (— l)^^^) only cause phase shifts, and do not 



alter the probabilities Pm^r 



Indeed, all functions 



f{x) lead to the same T2 relaxation rate in the inter- 
mediate 129 steps of the algorithm. This T2 relaxation 
rate should be compared to the overall T2 relaxation rate 
7/70 = 4 for the case that the n = 8 spins are in the same 
spin coherent state but interact individually with an ex- 
ternal reservoir. Note that a macroscopically separated 

linear superposition state ^(| -f + -- -+) + | )) 

has the same collective relaxation rate 7/70 = 4. 

One may infer that the Deutsch-Jozsa algorithm for 
this particular f{x) function is vulnerable to collective 
Ti relaxation processes, but is rather robust against col- 
lective T2 relaxation processes. 

Next, let us consider a different f{x) that takes value 
one when the last four qubits are of odd parity, and zero 
otherwise. In this case, the final state of the algorithm is 
IV'4) = 1+ + + + ). 

Figure IHl shows the evolution of P{j,m) of the n = 8 
qubit control register state. Even though the state \ip2) 
immediately after the first Hadamard transform is still 
subject to superradiant decay, the state |i/'3-3) imme- 
diately before the second Hadamard transform and the 
final state \1Jj4) have a large support on more stable sub- 
radiant states \j < n/2, m = —j).The decreased Ti relax- 
ation rate for this case, compared to the previous case, 
is obvious in Fig. [7j3. 



2. Improved implementation 

It is suggestive from the above examples that balanced 
functions f{x) which result in final states with m ~ 
have rather stable IV'3) and IV'4) states against collective 
Ti and T2 processes. Therefore, external longitudinal 
and transverse relaxation occurs primarily after the first 
Hadamard transformation. 

This problem can be resolved by starting the Deutsch- 
Jozsa algorithm with a state which has half up-spins and 
half down-spins. For instance, consider the initial state 

I h -I- -I- for the n = 8 case. This choice of the 

initial state does not change the fundamental structure of 
Deutsch-Jozsa algorithm at all; if the final state is found 
to be identical to the initial state, f{x) is a constant 
function. However, the modified initial state leads to a 
dramatic improvement in the Ti relaxation rate, as shown 
in Fig. [7p,D. 



B. Grover's data search algorithm 

1. Standard implementation 

One application of the Grover iterant in the database 
search algorithm consists of the following four steps, as 



shown in Fig. I^pl|. 

Stepl : Initialization of n qubit control register. 

I^i> - |0> = I > 

Step2 : Hadamard transformation. 

1^2) = m = -1= ^ |x) 

Step3 : Phase flip of the target state |t). 
t . |V3) = (/-2|T)(r|)i?|0) 

Step4 : Hadamard transform. 

Jl^ \yj,)=H{I-2\T){r\)H\0) 

Step5 : Phase flip of the |0) component. 

To . m = -{I- 2|0)(0|)i?(/ - 2\t){t\)H\Q) 

Defining the iterant Q as —loHIrH, the state jV's) can 
be rewritten as 

= -(/ - 2|0)(0|)ff2|0) + 2{T- 2\0){0\)H\T){r\H\0) 
= {l~4\Urfm + 2Ur,,H\r), (28) 

where Uro^ {t\H\0). 

The same rotation operator Q translates the state H\t) 

to 

QH\t) = -ToHTrHH\T) 

= {I-2\0){0\)H\t) 

= H|t)-2C/;J0). (29) 

It is clear from Eqs. (|28|) and l|29|) that the rotation 
operator Q preserves the two-dimensional vector space 
spanned by |0) and H\t). The operator Q rotates any 
linear superposition of |0) and H\t) by an angle 9 ~ 
2|f7T-Q| in this plane. The number of iterations that are 
required to transform an initial state |0) to a target state 
iJ|r) is therefore given by 



It 



(30) 



where |C/ro| — lv2^ is used. The target state is obtained 
with near unity probability by the Hadamard transfor- 
mation of the final state, H{H\t)) = \t). 

Each iterant can be split into the phase flip operation of 
the target state I^- and the inversion-about-the-average- 
probability-amplitude operation —HIqH We will 
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examine the system states after the inversion operation 
for the l^^ iterant. 

Fig. [TUl shows P{j, m) for n^8 after / = 1, 3, 5, 7, 10 
and 12 iterations, for the case of |t) = |255) = | + + + 
+ + + + +). Fig. ^] plots Q{9,ip) for the same states. 
As expected, both plots of P{j, m) and Q{9, (p) indicate 
a monotonic decrease in the probability of finding the 
initial state -ff|0) and a monotonic increase in the prob- 
ability of finding the target state. In this example, both 
the initial state -ff|0) and the target state |r) are in the 
completely symmetric (j — n/2) subspace, and so are all 
the intermediate states belong to this subspace. There- 
fore, one would expect that the implementation of the 
Grover algorithm for this target state is vulnerable to a 
collective longitudinal relaxation process. 

The solid line in the left panel of Fig. E| shows the 
normalized Ti relaxation rate at each of 50 selected evo- 
lution times, corresponding to the initial and final states, 
and the state after each gate in the 12 iterants. The 
Ti relaxation rate is highest immediately after the first 
Hadamard gate and slowest immediately after the second 
Hadamard gate. The oscillatory behavior per iteration 
continues to the midpoint (1 — 6 Grover iterations in this 
case). The high Ti relaxation rate corresponds to the 
state close to H\0) and the low Ti relaxation rate corre- 
sponds to the state close to |0). At the mid-point, the 
phase of this oscillation is flipped. The high Ti relaxation 
rate now corresponds to the state close to H\t = 255) and 
the low Ti relaxation rate corresponds to the state close 
to |255). 

In the left panel of Fig. ^1 the solid line indicates 
the normalized T2 relaxation rate at the same 50 com- 
putational steps. The relaxation rate is increased af- 
ter the first Hadamard gate since that state -ff|0) has 
a fairly large dispersion on the quantum number m, but 
is decreased after the second Hadamard gate, where the 
state is approximately ~ |0) — |iJ|255) (unnormalized) 
which has a much smaller dispersion on m. The T2 re- 
laxation rate is further increased immediately after the 
third Hadamard gate, where the state is approximately 
~ iJ |0) -I- 1 1 255). The T2 relaxation rate is high again im- 
mediately before the last Hadamard gate, where the state 
is approximately ^ -ff|255) which has the same disper- 
sion on m as the state H\0), and is low again immediately 
before the second last Hadamard gate, where the state is 
approximately ~ |-ff|0) -|- |255). 

Evidently, this implementation of the search algorithm 
is vulnerable to both longitudinal and transverse relax- 
ation processes. 

In Fig. m P{j,m) for n = 8 after I = 1,3,5,7,10 
and 12 iterations is plotted for a different instance of the 

search algorithm, where |t) = |15) = | |--|--l--l-)is 

the target state. As shown in the last panel , the target 
state |t) in this case has support largely in the subspaces 
where j € {2, 1, 0}. The state transforms from an unsta- 
ble superradiant state to a more stable subradiant state 
as the number of Grover iteration increases. The solid 



lines in the right panels of Figs. ^] and E| present the 
calculated Ti and T2 relaxation rates for this instance. 



2. Improved implementation 

The vulnerability of the Grover algorithm implemen- 
tation to collective longitudinal and transverse relax- 
ation processes stems from the initial spin-coherent state 
H\0). In order to circumvent this weakness, we can start 
with an alternative initial state such as I7 = 15) = 

I h + + +). If we replace a standard Grover 

iteration shown above with 

Q = -{T- 2|7)(7|)ij(/ - 2\t){t\)H , (31) 

this new rotation operator preserves the two dimensional 
vector space spanned by I7) and H\t). Ur^ — {t\H\^) is 
equal to l/\/2" or — 1-\/2", depending on I7) and |t). The 
plus and minus signs correspond to the opposite rota- 
tion directions of the initial vector by a Grover iteration. 
Here an absolute rotation angle 9 ~ 2.\Ut-^ \ = 2\{t\H^) \ is 
identical to the previous rotation angle 2Ura — 2(T|i?|0). 
Therefore, the ^ \/2" successive operations of Q rotates 
an initial state I7) to the final state H\t)^ from which we 
can find the target state |t) by a single Hadamard trans- 
formation. In this way we can avoid the unstable spin 
coherent state H\Q) from the implementation of Grover's 
algorithm simply by using any computational basis state 
with m = as an initial state. 

The improved Ti and T2 relaxation performance by 
this choice of an initial state is demonstrated by the 
dashed lines in Fig. 1121 and [T^ 



IV. CONCLUSION 

We have studied the collective decoherence proper- 
ties of n qubit quantum registers under implementation 
of selected quantum algorithms. A perturbative treat- 
ment of longitudinal relaxation implies that it is prefer- 
able to modify quantum algorithms to maximize the evo- 
lution time spent in the vicinity of subradiant states 
|j, m = — j), and avoid unstable superradiant states. 
Here, we have shown that slight modifications of the 
initial states of the Deutsch-Jozsa and quantum search 
algorithms are advantageous in weighting more heavily 
the subradiant states, if error correcting codes or more 
complicated logical encoding is prohibitive due to limited 
resources. Ideally, one might use the Ti transition rate 
as a parameter to optimize the gate decompositions and 
initial states used to implement a given algorithm. The 
optimal strategy whereby this can be done remains an 
open question. 

In our approach, the transition rates were calculated 
only at selected discrete evolution times in each algo- 
rithm. In a forthcoming paper, we will present a finer 
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grain analysis of the transition rates, where the algorithm 
is described in terms of single and two-qubit gates, and 
the state between each gate is considered in the over- 
all transition rates. In addition, an argument will be 
espoused as to how such an analysis can approximate 
a continuous average over the evolution generated by a 
time-dependent Hamiltonian consisting of pairwise inter- 
actions and external single-qubit pulses. 

The work of S.U. is partially supported by CREST, 
JST and that of C.P.M. and Y.Y is partially supported 
by NTT Basic Research Laboratories, SORST, JST and 
AFOSR under the contract of F49620-01-1-0556-P003-4. 



its inverse U ^ from right and left of the non-diagonal 
matrix, where 

X! ("^l,TO2,•••m^r|j,m, a). (A3) 

mi ,mjv 

Each column of this real unitary matrix provides the 
probability amplitude c„ij_...,„i^ for linear expansion of 
|j, TO, a) in terms of |mi, TO2,- ■ ■ , rrin)., i.e. 

\j,m,a)= ^ c„ii...„„|toi,TO2,- • -TOn). (A4) 



The coefficients Cn 



are the Clebsch-Gordan coeffi- 



APPENDIX A: CLEBSCH-GORDAN 
COEFFICIENTS 
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TABLE I: A diagonalizing matrix U defined by ljA3|l for the 
n = 3 spin case. The Clebsch-Gordan coefficients for {j, m) — 
(i' ~l) states represent the subradiant states. 



Since a symmetrized spin state \j, to, a) satisfies the 
eigenvalue relations in Eq. JSj) and the orthonormality re- 
lation, (j', to', a'\j, TO, a) = SjjiSmm'Saa' , the squared to- 
tal angular momentum is a diagonal matrix in this 
basis: 



{j',m',a\J \j,m,a) ^ j{j + l)SjfSjn,n'Sa 



(Al) 



By inserting two identity operators / = 
Smi- m„ l"^i'™2- • •TO„)(TOim2- • -TOnl bcforc and af- 
ter in the above equation, we obtain 

{j',m',a'\m[,m'2,--- ,m'^) 

{m[,m2,- ■ • ,to^|P|toi,TO2,- • • , m„) 
(TOi,m2,- • • ,TO„|j, TO, a) = j(j + l)djj'S (A2) 



This equation suggests that the non-diagonal matrix 
{m'l, TO2,- • • , TO^I J^|toi, TO2,- • • , toat) can be transformed 
into the diagonal matrix (j', m', a'\J'^\j, m, a) with eigen- 
values j{j -|- 1) by multiplying a unitary matrix U and 



cients |2ll |. 

According to the above argument, the symmetrized 
states \j, TO, a) are mathematically constructed by the 
following procedure. 

1. Divide computational basis states into the groups 
with different total longitudinal quantum number 
m = X]r=i"^i' -'^O'" instance, n — 3 spin states 
are split into four subspaces 

shown in the right column in Table 

2. Calculate in terms of computational basis states. 
Here, non-zero elements appear only in block diag- 
onal sub-matrices belonging to the same subspace 
designated by to. 

3. Find the real unitary matrix U which diagonalizes 

, where the diagonal elements of the diagonalized 
matrix are identically equal to the eigenvalue j ( j -I- 
!)■ 

4. The probability amplitudes of computational ba- 
sis states for constructing each symmetrized state 
(Clebsch-Gordan coefficients) are given by each col- 
umn of the real unitary matrix U. The example for 
n = 4 is listed in the left part of Table |I| 

The n-spin symmetrized states are grouped into dif- 
ferent subspaces by their eigenvalues j and m as shown 
in Table The first column with i = f represents 
a set of completely symmetric states, which is often 
referred to as angular momentum eigenstates. These 
states have no degeneracy, so the total number of the 
completely symmetric states is n -I- 1, corresponding to 
TO = f't ~ l''--7— f + The second column 

with j — f — 1 and m = ^ — l,---,—^ + 1 rep- 
resents d = nCi — nCo — {n — l)th fold degenerate 
states. The total number of eigenstates in this sub- 
space is (n — 1) X (n — 1) = (n — 1)^. The third col- 
umn with j = f - 2 and m= f-2, •••,-f-|-2 rep- 
resents d = „C2 — nCi — ^{n — 3)th fold degenerate 
states. The total number of eigenstates in this subspace 
is (n - 3) X f (n - 3) = f (n - 3)^ The final j manifold 
ends with either j — ^ (if n is an odd number) or j = 
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TABLE II: A distribution of symmetrized states \j, m, a) 



(if n is an even number). In_general the degeneracy of 
the states \j,m) is given by flj|, 



n!(2j + l) 



(t+J + l)!(t-j)! ' 
which is independent of m. 



(A5) 



APPENDIX B: COLLECTIVE LONGITUDINAL 
(Ti) RELAXATION RATE 

The collective spin-reservoir interaction Hamiltonian, 
which is responsible for a Ti relaxation process, is given 
by Eq. (Q). We will solve the Schrodinger equation, 



with an initial state 

|V(0)) = |j,m,a)®n,|0) . 



(Bl) 



(B2) 



We assume that all bosonic modes with a continu- 
ous spectrum are in ground states at i = 0. Ho = 
hojJ^ + hojialcii is a free Hamiltonian for the com- 
bined spin-boson system. A time dependent solution of 
the Schrodinger equation is 

I^W) =c™,o(t)e~*|j,m,a)|0) + 
^c,„_i,i.(i)e"^-(™'i)*-*-'*|j,m-l,a)|l,) , (B3) 

i 

where |0) =n,|0), and = |0)i • • • |0),_i|l), • • • . 
Substituting Eq. IjBSp into Eq. (jBip . we obtain 



«/i[c™,o(i)e~*b-,m,a)|0) 

i 

= fc4c™,o(i)e*"'"Vi(j + 1) - m{m " 1) 

i 

X \j, m-1, a)\h) + c™_i,i, 

X + 1) - m{m - m, a) |0)] . (B4) 



Multiplying {j,m,a\{0\ and {j,m,a\{li\ on both sides of 
Eq. HB4() . we have the coupled mode equations for Cmfi{t) 
and c,„_i,i,(t). 



Cm,o(t) 



Cm-i,ii(i) = -ikiV j{j + 1) - - 1) 

X e^(— ')*c„,o(t) , 

The initial condition for the above equations is 

Cm,o(0) = 1. 
c™_i,i,(0) = 0. 



(B5) 
(B6) 

(B7) 

We can integrate Eq. ljB6|l . taking initial condition Eq. 
(|B7|) into account, and substitute Cjn-i,ii{t) into Eq. 
(|B5|) to obtain the following integro-differential equation 

di'e^("-"')(*-*')c™,o(i')rfi' (B8) 



If we replace Cm,oit') in the righthand side of Eq. ljB8|) 
by Cm,o(0 replace the summation by an integral 
j p{ijji)duji with an energy density of states p{uJi), l|B8|l 
is reduced to 123 



(B9) 



where 



7 = + 1) - m(m - l)]p(c.) , (BIO) 

where we neglect the frequency shift. Since 70 — 
^k'^p{uj) is a longitudinal relaxation rate for a single 
spin, the collective relaxation rate is enhanced or sup- 
pressed by a factor of j(j + 1) — m(m — 1). 



APPENDIX C: COLLECTIVE TRANSVERSE (T2) 
RELAXATION RATE 



The collective spin-probe interaction Hamiltonian, by 
which a QND measurement of is realized, is given by 
Eq. (Q) . We assume the use of a single bosonic mode as a 
quantum probe and thus the summation over mode index 
i is suppressed. We will solve the Schrodinger equation 



ih-\i,it)) = 'Hi\m) 



(Cl) 



where a free Hamiltonian Tig is suppressed. The unitary 
time evolution operator in this case is reduced to 



U{t) = expl-igtJ'^a^a] 



(C2) 
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We assume the initial state for the spin and probe sys- 
tems is given by 

\i>{0)) =[cm|j, Ti, a) +Cm'\j',rn\a')]s 

«)^c„|n)p. (C3) 

n 

The state after the interaction is obtained by 

= ^[c„,e*9™"*|j,m,a), (C4) 

n 

+ c™,e^3'"'"*|j',m',a'),]®c™|n)p, 

where J^\j,m,a) and a^a|n) — n\n) are used. 

The reduced density operator 'pf'^'^'^ for the spin system 
corresponding to Eq. (|C4|) has the following non-zero 
components: 

p<-'^)^Tr,(|^)(^|) (C5) 
= [\cni\'^\j,m,a){j,m,a\ 
+ |c™f |j',m',a')(j',m',a'| 
+ Crac:^, \cn\'e-^^^">'\j, m, a){f, m', a'\ 

n 

+ Cm'C J2 |c„pe-^«(™'-™)"*|j', m', a'){j, m, a|] 

n 



Eq. IIC5|l suggests that the ofF-diagonal element has 
an oscillation frequency broadening Au = g\m — 
m'\{Any^^ where {AnY^^ = En |cn|'(n - (n))^ is 
the variance of the particle number. If we consider this 
frequency broading as a continuous variable, the off- 
diagonal elements c„iC*„/ and c„i'C^ relax with a rate 
given by 

r = cg\Tn-m'\{An^)i , (C6) 

where c is a constant determined by the distribution 
I Cm p. For a single isolated spin in a linear superposi- 
tion state where |m — m'| = 1, the transverse relaxation 
rate is Iq — cg{An'^)^ . Therefore, the collective trans- 
verse relaxation rate is either or suppressed by a factor 
of |m — m'l, compared to that of a single spin. 
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FIG. 1: Normalized matrix elements for upward (left) and downward (right) transitions as a function of j and m. 



z 




FIG. 2: Spin coherent states for arbitrary j manifolds. 




FIG. 3: Q{9,ip) representation of five states seen from the "north-pole" (top) and from the equatorial direction (bottom). 
Q — 1 appears as black, while Q = appears as white. The five states (unnormalized) from left to right are \ j = 4, m = 4) + 
\j = 4, m = -4), \j = 4,m = 3} + \j = 4,m= -3), \j = 4,m = 2} + \j = 4,m = -2), \j = 4,m = 1} + \j = 4,m = -1) and 
li = 4,m = 0>. 
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FIG. 5: Projections P{j, m) onto the symmetrized states for the Deutsch-Jozsa algorithm; n = 8 and f{x) is one iff the parity 
of X is odd. The plots, proceeding left-to-right in the first and second rows, correspond to the projections of \ipi) (top left) 
through I7/J4) (bottom right), indicated in Fig. |1| Note that the initial and final states for this example are \tpi) — |0) and 
IV'4) = |2"-1). 




FIG. 6: Q{6,ip) for the Deutsch-Jozsa algorithm; the plots correspond to the states described in Fig. |^ Note that in each 
graph, the unit sphere has been mapped onto the x-y plane; the transverse axes vary between < 9 < tv and — tt < ip < n. 
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FIG. 7: The normalized Ti relaxation rates at each of the 131 computational steps in the Deutsch-Jozsa algorithm. The traces 
correspond to: (A) Unmodified algorithm, f{x) = 1 if the parity of x is even. (B) Unmodified algorithm, f{x) = 1 if the parity 
of the last four bits of x is even. (C) Improved algorithm, f{x) as in A. (D) Improved algorithm, f{x) as in B. 




FIG. 8: P{j,m) for states \tpi) through \tp4) in Doutscli-Jozsa algorithm, n — 8, f{x) equalling one if the parity of the last four 
qubits is odd. Note that the initial and final state are \ip-i) = |0) and |'i/'4) = |15). 
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FIG. 9: Circuit topology for Grover's database search algorithm 




FIG. 10: P{j, m) of the control register states in the n = 8 database search algorithm. The plots correspond to the states after 
Grover iteration i G {1, 3, 5, 7, 10, 12}; ^ = 1 is top left and ^ = 12 is bottom-right. The initial and final states are = |0) 
and IV'4) = |2" - 1). 




FIG. 11: Q{6,ip) of the control register states in the n = 8 database search algorithm. The plots correspond to the states after 
Grover iteration £ € {1, 3, 5, 7, 10, 12}; £ = 1 is top left and ^ = 12 is bottom-right. 
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FIG. 12: The normalized Ti relaxation rates at each of the 50 computational steps for the n = 8 database search example; 
the left figure applies to the target state |r — 255) and the right figure is for |r — 15). The first Hadamard gate is followed by 
twelve Grover iterations consisting of four steps of —HIoHIt- The solid and dashed lines represent results for the standard 
implementation using an initial state = |0) and the improved implementation using an initial state |^i) = |15) = 

I 1- + + +), respectively. 





FIG. 14: P{j, m) of the control register states for the second n = 8 example of the database search algorithm. The initial and 
final states are |?/>i) = |0) and \^l)4,) = |15). 



